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Turbulence in the Intergalactic Medium 
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ABSTRACT 

We study supernova-driven galactic outflows as a mechanism for injecting turbulence 
in the intergalactic medium (IGM) far from galaxies. To this aim we follow the evolu- 
tion of a 10 13 M Q galaxy along its merger tree, with carefully calibrated prescriptions 
for star formation and wind efficiencies. At z ~ 3 the majority of the bubbles around 
galaxies are old (ages > 1 Gyr), i.e. they contain metals expelled by their progen- 
itors at earlier times; their filling factor increases with time reaching about 10% at 
z < 2. The energy deposited by these expanding shocks in the IGM is predominantly 
in kinetic form (mean energy density of 1 [ieV cm" 3 , about 2 — 3x the thermal one), 
which is rapidly converted in disordered motions by instabilities, finally resulting in 
a fully developed turbulent spectrum whose evolution is followed through a spectral 
transfer function approach. The derived mean IGM turbulent Doppler parameter, bt, 
peaks at z w 1 at about 1.5 km s" 1 with maximum b t — 25 km s" 1 . The shape of the 
bt distribution does not significantly evolve with redshift but undergoes a continuous 
shift towards lower b t values with time, as a result of bubble aging. We find also a 
clear trend of decreasing bt with Nm and a more complex dependence on R s resulting 
from the age spread of the bubbles. We have attempted a preliminary comparison with 
the data, hampered by the scarcity of the latter and by the challenge provided by the 
subtraction of peculiar and thermal motions. Finally we comment on the implications 
of turbulence for various cosmological studies. 

Key words: cosmology: theory - galaxies: intergalactic medium - galaxies: evolution 
- quasars: absorption lines. 



1 INTRODUCTION 

The intergalactic medium (IGM) is a pervasive, diffuse cosmic 
baryonic component out of which galaxies form by accretion. In 
turn, the IGM is replenished with gas and heavy elements carried 
by galactic outflows. Such dynamical processes hence encode a 
record of the complex galaxy-IGM interplay. 

Observations of metal absorption lines (e.g. C III , C IV, 
Si III, S IV and O VI) in quasar spectra are a direct probe of 
this evolution. They show that regions of enhanced IGM density 



significant increase of the C IV abundance tow ards lower red- 
shifts z ,§ 2. This result has been also confirmed bv lCooksev et all 
j201Ch for z < 1. 

Since metals are created in stars inside galaxies, their 
transport into these distant regions must rely on some yet 
unknown mechanisms. The most likely one i s supernova- 



driven galactic outflows llMac Low fc Ferraralll999l;l Ferrara et al 
| 200C| [ hereafter FPS1; lAguirre et al.ll200lt lOppenheimer fc Dave 
2006; ICen fe Chisaril l2010h . although different alternatives 



far from large galaxies are polluted with non-neglif 
of heavv elements (e.g. ISongaila fc Cowiej Il99£ 


;ible amounts 
; iDave et alj 


199SI; ISchave et alj 


20001; iPettini et alj 20031; ISchave et al J 120031; 


Aguirre et alj2004 


Aracil et alj 2004; Simcoc et al. 2004). In ad- 



dition, hya forest clouds with neutral hydrogen column density 
log N(HI) > 14 are metal enriched to Z & 10~ 3 - 10~ 2 Z B 
JCowie et alj 1 19951 ; ISongaila fc Cowid Il996l ; lEUison et a"l]|200Ch ; 
finally, the carbon and silicon abundances of such systems remain 
roughly constant throughout the redshift range 1.5 < 2 < 5.5 
JSongailall200ll ; IPettini et alj|2003l) . although the first signs of a 
slight decrease above z = 6 have been now tentatively reported 
jRvan- Weber et alj [200$) . More recently, at varian ce with the 
constant behaviour picture, iD'Odorico et al" . (2010) reported a 



E-mail: evoli@sissa.it 



have been proposed, i .e. dust-sput tering Aguirre et alj l200ll : 
iBianchi fc Ferrarall2005l) or merging jGnediri 1 19981) . This hypoth- 
esis is supported by a series of observations of galactic winds 
(e.g. iHeckmanl 12001 ; IPettini et alj|200ll; lAdelberger et alj |2003|; 
Shap lev et al.l2003l:lErb et al.l2 006) from which one can conclude 



that (on some spatial scale) galactic wind velocities range from 
hundreds to thousands of km / s and the mass loss rates are com- 
parable to star formation rates (SFRs). 

A complementary piece of information comes from the well- 
established fact th at the inters t ellar medium (ISM) in galactic 
disks is turbulent jLarsonlll98ll ; IScalolll987l ; iDickev fc Lockmanl 
1 199d : lElmegreen fc Scaloll2004l ) . In the Milky Way (MW), the ob- 
served vertical distribution of gas, and especially the cold neutral 
HI component, cannot be supported by thermal pressure alone 
jLockman fc Gehmanlll99ll ; lFerrarall993l) and an additional pres- 
sure contribution must be advocated. Furthermore, pulsar scin- 
tillation experiments have convincingly demonstrated that elec- 
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tron density fluctuations in the ionized ISM are characterized 
by a Kolmogorov spectrum on scales r anging from AUs to kpc 
llArmstrong et aL I ll995l ; lHan et alJkoCMh . The Kolmogorov spec- 
trum is the spectrum of the density fluctuations which arise from 
turbulent motions in a compressible gas with constant dynami- 
cal energy exchange between irregularities on different scales and, 
under very general conditions, it assumes the form of a power-law 
(~ A; -11 / 3 where k is the wave-number associated to a given scale 
L) within the inertial range. In other spiral galaxies the vertical 
velocity dispersions (derived essentially from HI observations) are 
observed to vary radially from 12 — 15 km/s in the central parts 
to 4 — 6 km/s in the outer regions and, in most cases, these val- 
ues exceed th ose expected by the thermal broadening of the HI 
emission line toib et al.ll2006f) . 

Often invoked to provide the required additional pressure 
are turbulent mot ions of the interstellar gas, magnetic fields, 
and cosmic rays |Bp ulares fc Cox Il990|; lLockman fc Gehmanl 

Turbulence is likely 



Il99ll; iDib et all I2OO6I ; Ijoun" et al 



to be to most relevant pressure contribution; its driving may 
be due to different sources (star formation, stellar outflows, in- 
stabilities) acting on specific characteristic scales and injecting 
different amounts of kinetic energy into the medium; among 
these, supernova (SN) explosions are thought to be the dom- 
inant agent, at least in the star forming regions of galaxies 
llMcKee fc Ostrikerll977l;lNorman fc Ferrar jl996tlMacLowl2004 
[de*^vT!leT^rFireitscIrwer^?l2007l) . If both turbulence and metal 
enrichment are predominantly driven by the same physical phe- 
nomenon, i.e. kinetic energy deposition by supernova shocks, then 
a natural expectation is that IGM gas in enriched regions should 
show signatures of a turbulent regime. 

A direct way to measure turbulence in the IGM is to look 
for velocity differences on the smallest spa t ial sc ales accessible 
to observations; to this scope, iRauch et ail j200lT) m a pioneer- 
ing experiment used lensed quasars. They observed adjacent CIV 
profiles in paired lines of sights, and find velocity differences of 
Sv ~ 5 km s _1 on scales of 300 pc at redshift ~ 2.7. They apply 
the Kolmogorov steady-state assumption to calculate an energy 
transfer rate e fa 10 -3 cm 2 sec -3 and a turbulent dissipation 
timescale as short as 10 s yrs, thus implying a recurrent stir- 
ring of the gas by some sort of galactic-scale feedback. 

Another hint of a significant turbulent contribution to the 
IGM kinetic budget comes from the fact that the median Dopplcr 
parameters measured in the Lya forest are significantly larger 
than those predicted by cosmological simulations. Again, this im- 
plies that some energy in non-thermal form must be injected in 
the gas to explain the observed line broadening l lMeiksin et al.l 
l200l];lMeiksiii(l2009l') . 

Quite surprisingly, the properties of intergalactic turbu- 
lence have received so far relatively little attention in spite 
of the insights that its detailed underst anding might provide 
on the galaxy-IGM interplay ( FPS; lOppenheimer fc Pavel 20091 ; 
IScannapieco fc Fj riiggcn 201jJ|. [Oppenheh*n*er**fc*rJaveri200£ ) try- 



ing to model low-z OVI absorbers concluded that their proper- 
ties fit observables if and only if sub-resolution turbulence (in 
practice a density dependent turbulent fe-parameter) is added at 
a level which increases with OVI absorber strength (and hence 
presumably closer to high density regions hosting galaxies). In 
addition, stronger absorbers arise from more recent outflows. It is 
important to note, that turbulence, if present, can also profoundly 
modify the enrichment patterns through diffusion processes, as 
pointed out already by FPS. Additional important information 
might come from investigations of structure formation, during 
which the released gravitational energy is transferred to the IGM 
in different forms (e.g. gas entropy, instabilities, cosmic ray ac- 
celeration); in particular, turbulence is induced by the vorticity 
cascade originating at cosmological shocks. Th is has been inves - 
tigat ed throughout c osmological simulations by I Ryu et al. I fcoosl ) 
and lZhu et al] l|2010h respectively, who derived the average mag- 



netic field strength and turbulent pressure in the overdense IGM 
regions outside clusters/groups. 

In this work, we try to fill the aforementioned gap by investi- 
gating in detail the supernova-driven outflow scenario as a mech- 
anism to pump and sustain turbulence in the IGM surrounding 
high-redshift galaxies. We investigate the evolution of winds via 
a semi-analytic approach following both galaxy evolution along 
its hierarchical growth and the expansion of supernova driven 
superbubbles as they escape from the halo potential well. With 
such results in hand, we model turbulence evolution within wind 
shells using a novel approach based on the spectral transfer equa- 
tion. This approach allows to derive the turbulent energy density 
deposited by galactic winds in the IGM, along with its spectral 
features and dissipation time scales, and to predict the corre- 
sponding thermal/kinetic evolution of the IGM. 



2 THE MODEL 

In this Section we describe how we determine the evolution of the 
turbulent energy of the deposited by galactic outflows into the 
IGM. First, we introduce the galaxy formation model, allowing a 
precise description of both the SF history and the SN feedback 
along the galaxy merger tree (MT). These results are then used 
as inputs to model the evolution of pressure-supported bubbles 
through the solution of dynamical equations and the turbulence 
energy spectrum evolution within them. 



2.1 Cosmic star formation history 

The semi-analytic model we introduce in this Section to describe 
the evolut ion of a galaxy along i ts MT is similar to the one-zone 
model by Salvadori ct al. (2007) with the improvements intro- 
duced by ISalvadori et al.l J200ST) . Since in this work we are not 
interested in following the chemical enrichment in detail, we use 
a simplified approach described in the following. 

In hierarchical models of structure formation, such as 
ACDM, the formation of a DM halo thro ugh accretion and re - 
peated mergers can be described by a MT jLacev &: Colelll993h . 
MTs, which list the progenitors of a given halo at different red- 
shifts and describing how and when these merge together, con- 
tain essentially all the necessary information about the dark mat- 
ter content of halos to build realistic model of galaxy formation 
as demonstrated by a lon g-standing practice in the field (e.g. 
iKauffmann fc Whitej[l993l; ISomerville fc KolattHl999T ). 

The MTs can be extracted from N-body simulations 
llSpringel et alj l2005t) or generated by Monte-Carlo algorithms 
llSheth fc Le mson 1999, e.g). In this work we compute the mass 
growth histories of DM halos using the second method using 
a pu blic co dttJ based on the extende d Press-Schecter formal- 
ism JCole et aij|2000l) , as improved by Parkinson ct al. 
to which we defer the reader in numerical details. As an input, 
we use a cosmological parameters and power-spectrum. 

As an output, we obtain many realizations of the MT for a given 
halo mass Mo at a certain final redshift Zf, each realization lists 
all progenitors of Mq at different redshifts, following the related 
merging histories down to a mass of M IeB . 

In our model all the proto-galaxies virializes with a gas-to- 
dark matter ratio of 1:5. The neutral gas in these mini-halos 
cannot cool via atomic hydrogen and relies on the presence of 



1 |http : //s tar-www . d ur ■ ac ■ uk/~] cole/merger_trees 

2 Throughout all the paper we adopt a flat ACDM cosmological 
model with h = 0.72, U b = 0.044, U m = 0.26, n s = 0.96 and 
erg = 0.8 consistent w ith the five-year WMAP parameter analysis 
llDunklev et al.ll2009h . 
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molecular hydrogen, H2, to cool and collapse, ultimately form- 
ing stars. Since we assume that feedback effects rapidly suppress 
star formation in the first minihalos and that only Lya cooling 
halos (T v ; r > 10 4 K) contribute to the SF history of the galaxy, 
we choose M TeB = M±(z = 20) where 4/4(2) is the mass corre- 
sponding to a halo with virial temperature of 10 4 K at redshift 
2, given 



M 4 (z) « 10 8 M Q 



10 



3/2 



(1) 



,1 + 2/ 

(for an exact expression see lBarkana & Locb (2001)) and 2 = 20 
is the starting redshift of our simulation. Note that the value 
of Mres sets the limit between progenitors and mass accretion, 
since the galaxies above M re s experience star formation and stel- 
lar feedback, ultimately changing their gas content, whereas ob- 
jects below this threshold retain their original gas content fraction 
during their evolution, which is finally inherited by another halo 
in the next hierarchy level. 

Star formation in gas clouds occurs on a free-fall timescalc 
iff = (3tt/32Gp) 1 ^ 2 where G is the gravitational constant and 
p is the average mass (dark+baryonic) densi ty inside the halo 
assumed to have a virial radius (r v ; r ) as given in lBrvan fe Normanl 
(1998). The star formation efficiency of a galaxy is then modeled 
as a fraction, e, , of the free-fall time: 



M e {z) 
' tg{z) 



(2) 



where e„, represents a free parameter of the model and M g repre- 
sents the gas mass in the halo which has not yet been converted 
into stars. 

According to the standard scenario for galaxy formation the 
gas inside the galaxy is depleted by various feedback processes, 
the most important being gas-loss driven by SN explosions. In 
fact, SN explosions may power a wind which, if sufficiently ener- 
getic, may overcome the gravitational pull of the host halo lead- 
ing to expulsion of gas and metals into the surrounding IGM. To 
model this process, we compare the kinetic energy injected by 
SN-driven winds with the minimum k inetic energy of a mass M w 
to escape the galactic potenti al well jLarsonlll974 ; iDekel fc Silk! 
ll986l : lFerrara fc Tolstovll2000l ). The mass of gas ejected from the 
galaxy is computed from the equation: 



where 



Est 



M w vl = E S n 



vN S n(Esn) 



(3) 



(4) 



is the kinetic energy injected by SNe and v\ = lE^jM is the 
escape velocity of the gas f rom a halo with mass and binding 
energy Ef, given bv lBarkana fc Loebl (120011) : 



1 GM 



2 R v 



M h 



5/3 



1 + 2 



5.45xlO j;i —5 2 — — ft" 1 erg. 

(5) 

In Eq.[4] e w is a free parameter which controls the conversion 
efficiency of SN explosion energy in kinetic form, N§n is the 
number of SN, and (Esn) = 1-2 X 10 51 erg is the average SN 
explosion energy (not in neutrinos). Differentiating with respect 
to time Eq.[3]we find that the gas ejection rate is proportional to 
the SN explosion rate: 



2ew(E S]v ) 



N st 



(6) 



In our model we only consider Type II SNe (pair-instability super- 
novae are neglect ed as PopIII stars do n ot contribute significantly 
to the total SFR (Tornatorc ct al. 200?])); hence we calculate iVgrj 
by integrating the adopted IMF (sec below) in the canonical su- 
pernova range 8-100 Mq. For any star forming halo of the galaxy 
hierarchy, we therefore solve the following system of differential 



equations: 
M,(z) = 

M g (2) = 



[1 - R(z)] i>(z) 

M a (z) - [1 - R{z)] V(2) - Af w f» 



(7) 
(8) 



The first equation simply defines the SFR. In our model, we 
have assumed the Instantaneous Recycling Approximation (IRA 
lTinslevll980T) . according to which stars are divided in two classes: 
those which live forever, if their lifetime is longer than the time 
since their formation and those which die instantaneously, even- 
tually leaving a remnant. The transition mass between the two 
possible evolutions, or turn-off mass (rrtto), has been computed at 
any considered redshift. Under IRA, the returned fraction, that 
is the stellar mass fraction returned to the gas through winds and 
SN explosions is: 



R(z) 



fm° (z) I m _ Wm(m)}<j>{m)dm 



r lOOM 

Jo.lMg rn4>(m)dra 



(9) 



where <f>(m) is the IMF of the newborn PopII/I stars and it is 
assumed to have the form: 

m cu t\ 



, s dN -l+x 
(m) = <x m ^ 

dm 



cxp ( 



111 



(10) 



with x = — 1.35, TOcut = 0.35 Mq and m in the range 
[0.1,100] M llLarsonl I 19981 ). The quantity w m (m) represents 
the mass of the stellar remnant left by a star of mass m 
which explodes as SN. We have u sed the grid of models by 
(Ivan den Hoek fc Groe ncwcgcn 1997) for intermedia te mass stars 
(0.9 < m < 8 Mm) and dWooslev fc Weaver! 1 19951) for massive 
stars (8 < m < 40 Mq). With the adopted IMF the corrispondent 
SN energy per unit stellar mass formed is 1.36 X 10 49 ergs/M . 

The second equation (Eq. 8) describes the mass variation 
of cold gas: the latter increases due to accretion (M a (z)) and 
it decreases both because of astration and mass loss due to SN 
winds. To model gas accretion we assume that if a new halo vit- 
alizes out of the IGM gas it forms with a cosmological gas fraction 
fb = r2b/f2 m = 1/6, whereas if it results form the merging of two 
already existing halos its gas content is simply the sum of the gas 
mass of the progenitors. Similarly, during the merger, the stellar 
mass of the new galaxy is assumed to be the sum of the stellar 
masses of the progenitors. The model free parameters (e w and 
e*) are fixed to match the global properties of the MW as we will 
discuss in more details in the next Section. 



2.2 Bubble evolution 

In our simulation galactic outflows are treated as pressure-driven 
bubbles of hot gas emerging from star-forming galaxies. They 
expand working against IGM pressure, and are driven by the en- 
ergy injected by multi-SN explosions. Most of the swept-up mass, 
both in the early adiabatic and in the following radiative phases, 
is concentrated in a dense shell bounding the hot over-pressurized 
interior. 

For this reason, galactic bubbles are canonically studied 
by using the thin-shell approximation in which the shell expan- 
sion is driven by the internal energy (Ej,) of the hot bubble gas 
dOstriker fc McKeeill988l: iMadau et aU200ll; iBertone et"al]|2005l: 
ISamui et al]|200Sl) . 

In this work we follow the modelization already presented 
in Mad au et al. I ll200ll ) (see their Sec. 2). However, in our model 
we derive the SFR evolution from a constrained galactic model 
and from the SFR we derive the mechanical luminosity, defined 
as L(t) = dE/dt where E is the energy produced by the total 
contribution of Ngfj SNe with energy {-Eg at) and an efficiency 

In the following, we remind the reader of the relevant equa- 
tions used to follow the evolution of the interesting bubble prop- 
erties, namely the bubble radius R s , the shock velocity V B and 
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the internal bubble temperature T b : 
d 



dt 



(VspRs) = 4nRi(P b - P) - 



GM(Rs) 
R% 



dE b _ 
dt 
dT b 
dt 



L(t) - A-nR 2 s P b R s - V s n 2 H!b A(f b ) 



23 Ci kT, 



,9/2 



R S P b 10 C 2 RjP b 



(11) 
(12) 

(13) 



where G is the gravitational constant; M(R) is the galactic DM 
mass assumed to follow a NFW profile within a radius R; p (P) 
is the density (pressure) of the ambient medium; V s = (47r/3)-?? 3 
is the volume enclosed by the shell; P b = E b /2-wR? s is the pres- 
sure of the bubble gas assumed to be isothermal and with adi- 
abatic index 5/3; n 2 H () A(f ) is the cooling rate per unit volume 
of the hot bubble gas (whose average hydrogen density and tem- 
perature are ftjf b and f, respectively); C\ = l6Trfim p r)/25k and 
C2 = (125/39)7r/im p (where rj = 6x 10 — 7 (c.g.s. units) is the clas- 
sical Spitzer thermal conduction coefficient, assuming a Coulomb 
logarithm equal to 30, and fi = 0.59 is the mean molecular weight 
of a primordial ionized H/He mixture). 

Finally, the cooling function (represented by A(T) in Eq. ll2jl 
depends on the hot bubble density, temperature and metallic- 
ity (Z). In principle, one could compute an average metallicity 
of the bubble, but for simplicity we assume for this quantity 
a constant value Z = O.IZq, noting that the final results are 
very slightly dependent on this choice within a reasonable range. 
The cooling rates due to ga s radiative processes are taken from 
ISutherland fc Dopital dl993h ; the other relev ant cooling agent is 
inver se Compton cooling off CMB photons lllkeuchi fc Ostrikerl 
1986), which is dominant at higher redshifts. 

The first two equations (Eq, llll and ll2H derive from momen- 
tum and energy conservation. In particular, the right-hand side 
of Eq. 1111 represents the momentum gained by the shell from the 
SN-shocked wind which expands against external pressure and 
gravitational attraction; the right-hand side of Eq. 1121 describes 
the mechanical energy input, the work done against the shell, and 
the energy losses due to radiation. The third equation (Eq. II 3D is 
obtained by equating the rate at which the gas is injected from 
the shell into the cavity with the conductive evaporation rate. 
Note that, as the bubbles sizes are always much smaller that the 
horizon scale, the cosmological terms in the above equations can 
be safely neglected. 

Another important quantity to determine the fate of the 
wind, is the total mass in the shell, M s , since the energy required 
to accelerate it increases with its mass. The shell grow rate M s 
is wind mass deposition rate (Eq. [3]l and the rate at which gas is 
swept- up from the galactic environment (for R s ^ rvir) or from 
the IGM {R s > Rvir)- 



- 4-irpR- 



,dR s 



dt 



(14) 



dA/ s 
dt 

The ambient gas density is taken to be equal to the halo 
gas density distribution (assumed to have an isothermal profile) 
within R v i r , and to the IGM background density outside the virial 
radius. To determine the pressure of the ambient medium we must 
further assume the ambient gas behave as an ideal isothermal gas 
which , inside halos, is at the virial temperature iNavarro et al.l 
(1997) and outside R vlr , the IGM, is taken to be photo- heated 
by the UV-background radiation to temperature calculated as in 
e.g. IChoudhurv fc Fer rara (2006;). Note that, even prior to reion- 
ization, the IGM into which the bubble expand will be locally 
ionized and heated to roughly the same temperature by the SN 
progenitor stars. 



3 This result has been derived by substituting the expression for 
the internal energy E t ^ = 3/2Nk b T into the ideal gas law P = 



However, we notice here that the bubble reaches very rapidly 
the virial radius, so that these bubbles are practically expanding 
in the IGM f o r mos t of their time. To see this, we consider as in 
iMadau et al. I fcooil i a typical objects virializing at high-rcdshift 
(z = 9) with an halo mass M = 10 s Mq. At these epochs, the 
dark matter halo of a subgalactic system will be characterized by 
a virial radius 

-1 



R vil = 0.76M„ 1/3 h 



1 + z 
10 



kpc 



(15) 



The Sedov solution (which is a good approximation of the equa- 
tion systems when the ambient gas pressure, gravity and cooling 
can be neglected) predicts that the shell radius evolves according 
to: R a = (125/1547T) 1 / 5 (Lt 3 /p) 1/5 ; assuming a constant lumi- 
nosity of 10 38 erg -1 , we find that the time taken to reach -R v i r is 
~ 10 7 yr. This time happens to be smaller than the Hubble time 
at that rcdshift. 

The velocity of the shock front decreases during the time be- 
cause of adiabatic expansion and cooling of the bubble, it might 
happens the shock velocity eq uals the inter g alacti c gas sound 
speed, in this case, similarly to lBertone et all (120051 ). we assume 
the dynamics of the wind joins the Hubble flow and no more mass 
is accreted. 

Since the Eqs. 1111 1121 and 1131 requires an initial condition 
which cannot be (no bubble) to be consistently solved, we used 
the Sedov solution for a very short time with respect to tji(z = 
20) to obtain the primordial bubble to follow the evolution. The 
evolutionary equations (Eqs. 1111 1121 and I13H can be integrated 
numerically to follow the evolution of the shell along with the 
thermodynamic properties of the bubble for a given halo. 

A final point concerns the treatment of bubbles when two 
halos merge along the hierarchical merger tree. In that case, we 
impose mass conservation, i.e. the mass the new shell is the sum 
of the masses and M 2 of the two progenitor bubbles, that is 
m/ = Ml + M s 2 . Similarly, the final volume is taken to be the 
equal to sum of the volumes of the two single cavities V/ = V-[ + 
V 2 , and by the adopted spherical symmetry, we update the shock 



/ 



[(Rl) 3 + (R 2 S ) 



3lV3 



Finally, the shock velocity 
/ _ 



radius as R 

is given by the momentum conservation, which states that v\ 
M^vl+M 2 v 2 where MJ = Mi/(M} + M 2 ), and the new internal 
energy is the sum of the progenitor bubbles internal energies, that 
is = E^ + E 2 . The gas temperature in the merged bubble is 
computed by assuming the final thermal pressure is the sum of 



the thermal pressures of the merging bubble: T, 



f 



where 



Pb 



pI/(pI+pD- 



-Pb j 



-Pb 



T 2 



2.3 Turbulence evolution 

In order to model the evolution of the turbulence developed in 
the expanding shells we adopt in the following an approach based 
on th e spectral transfer equation derived by iNorman fc Ferrara! 
l|l996T ) based on the hydrodynamic Kovasznay approximation. In 
fact, our aim is to derive averaged properties of the turbulent 
spectrum. 

It is useful to introduce a spectral representation of the ve- 
locity field in the fluid. Then the energy density in eddies with 
wave-number between fc and fc + dk is given by p£(k, t)dk, where 
p is the density of the fluid and £(k,t) is turbulence spectrum. 

The time evolution of this quantity, under the assumption of 
local mode interaction in every region of the spectrum, is given 
by the following equation: 
d 



a— [f(fc,t) 3 / 2 fc 5 / 2 ] - 2uk 2 £(k,t) +S(k,t) (16) 



nRT/V. 



where v is the kinematic viscosity and S(k,t) is the source func- 
tion and a is a dimensionless constant that should be fixed ex- 
perimentally. 

In order to describe the evolution of the turbulent spectrum, 
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we have finally to provide the appropriate boundary and initial 
conditions: 



£(0, t) = £(oo, t) = , < t < oo , 
£(fc,0) = £ (fc) , ^ k < oo . 



(17) 
(18) 



A caveat must be made here. The assumption underly- 
ing the above description is one of incompressible, homoge- 
neous, isotropic turbulence. This assumption clearly fails in the 
highly compressive IGM. However, the bulk of our knowledge on 
turbulence comes from terrestrial laboratory experiments, and 
most of them deal wi th liquids; in add i tion, s imulations of com- 
pressible turbulence jLazarian fc Choi [20(53 ; iKowal fc LazariarJ 
1 20071 ; ISchmidt et al,l2009l) have shown that nonlinear interactions 
rapidly transfer most of the energy to non-compressible modes, 
in final one can use those results at least as a reasonable guide 
when discussing compressible turbulence. 

To model S(k,t) we assume that turbulent motions in the 
outflow arc induced by interacting blast waves. That is, we main- 
tain that observed turbulent motions in the shell are ultimately 
derived from the kinetic energy of SN induced shock waves which 
can act as source function for the turbulent cascade with an ef- 
ficiency of order unity. If this is the case, the re are some general 
constr aints on the form of the source function. iBvkov fc Toptvginl 
(1987) have studied in detail the shock-induced turbulence phe- 
nomenon, and conc l ude t hat the turbulent spectrum has, for the 
iMcKee , 

sions, a k~ 2 dependence in the short-wavelength regime, while 



Ostrikerj j 19771) model for the SN shock-wave expan- 

tigth regime, while 
Thus, the simplest 



for long wavelengths it is proportional to fc 2 
source function retaining this behavior is 

S(k,t) = S (t)f(k/k ) (19) 

where f(x) = x 2 /(l+z 4 ), kg is the wave-number corresponding to 
the characteristic length at which the turbulence is injected (and 
assumed to be kg ~ l/R s ) and So (t) is a normalization factor ob- 
tained by equating the fe-intcgral of th e source function to the ki - 
netic energy rate (-Ekin)- Finally, as in Brandenburg et al.l (120071 ) 
we use a constant kinematic viscosity of v = 5 X 10 24 cm 2 s _1 ; 
this is equivalent to scale the dynamic viscosity coefficient with 
density. 

The spectrum evolution is followed for each halo of the MT; 
again, we need to specify how we treat merger events. To conserve 
turbulent energy associated to any MT node we assume 



which entails 



p£(k) = pi£l(fc)+/9i£ 2 (fe), 



£(*) = ^£ l( *) + ^£ 2 (*). 



MF? A "^" J ' MR 3 sl 
The turbulent pressure is finally given by 



Pt(t) 



f' CO 

> / dk£(k, 
Jo 



t). 



(20) 
(21) 

(22) 



3 RESULTS 

In this Section we present the results of our model. After a de- 
scription of the model calibration, useful to fix the model free 
parameters, we give our prediction for the bubble evolution along 
a MT extending down to z = along with the IGM turbulence 
properties. 



3.1 Model calibration 

Our galaxy evolution model includes several relatively poorly 
known (albeit important) physical processes that need to be em- 
pirically calibrated. To this aim we have used the observed prop- 
erties of the MW as a benchmark to fix the best values for the 
two model free parameters: e*, the efficiency of star formation 
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Figure 1. The 68%, 95%, 99% confidence level regions of our 
galaxy model are represented in the (e w — e„) plane. 



(Eq. [2]l and e w , the SN feedback efficiency (Eq. [3j. This has been 
accomplished by producing the MT of a dar k matter halo w ith 
mass consistent with the MW one (10 12 M^. IXue et alj|2008l) at 
redshift z = 0. We have then compared the resulting properties of 
the synthetic galaxy with the following ones deduced from MW 
observations: 

• Stellar mass. Contributions to the stellar mass come from 
the disk Mf ak ~ (4 - 6) x 10 1(1 M Q , the bulge M* u ' 9e ~ 
(0.4 - 1) X 1 10 M(7), and the halo Mj°- l ° ~ (0.2 - 1) X l O w Mp, 
components (Dclmcn & Binncy 1998: Brown et al. 2005), yield- 
ing a total stellar mass of M, ~ 6 X 1O 1O M . 

• Gas-to-stellar mass ratio, M gaB /M t 



0.13. The mass of 



gas has been derived using the observed mass of HI and HII re- 
ions of the Galaxy, M aas = M HI + Mrii ~ (6 + 2) X 10 9 M Q 
Stahler fc Pallall2005h . 
• Current star formation rate. Murray fc Rahman! l|2010T l use 

the total free-free emission in the WMAP foreground map as a 

probe of the massive star population and derive a global SFR of 

1.3 M yr" 1 . 

For all these observables we assume a relative uncertainty 
of 20%. Fig. [T] shows the \ 2 confidence levels of our model with 
respect to these observables as function of e„ and e w . The distri- 
bution presents a clear minimum at (e„, e w ) = (0.5, 0.012); how- 
ever, a degeneracy exists with a relatively wide range of choices 
that can reproduce with sufficient accuracy the global MW prop- 
erties. Thus, the best-fit model implies a star formation timescale 
which is only a factor 2 longer than the free-fall time and a 
relatively inefficient feedback in order not to expel from very first 
proto-galaxies too much gas which will then be crucial to fuel the 
subsequent SF in larger objects. 

From the same run we can obtain the redshift evolution of 
the SFR density along the MW merger tree for our best-fit model, 
after averaging over 100 different realization of the MT0. For com- 
parison sake only, we overplot the cosmic SFR on top of the MW 
SFR history. Perhaps not coincidentally, the two evolutions match 
each other quite well indicating that the Galaxy is a not too biased 



4 The effective comoving volume of the MW is taken equal to 
30 Mpc 3 , and corresponds the size of a 10 12 Mq linear fluctuation. 
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Figure 2. Redshift evolution of the SFR density along the MW 
merger tree: (top) for our best-fit model (e* = 0.5, e w = 0.012), 
(bottom) for the no-feedback model (e, = 0.5, e w = 0). The 
curves are obtained after averaging over 100 realizations of the 
MT; shaded areas denote ±1<j dispersion regions around the 
mean. Point s represent the low-r cdshift measurements of the cos- 
mic SFR bv lHopkins fc Beacoml J200d) . 



tracer of average cosmic conditions. Finally, in the same figure we 
show the SFR density evolution for a model with e w = 0; you can 
see that in this case the model is far from the data, which means 
that, even if the efficiency happens to be small, the feedback by 
galactic winds is crucial to obtain a consistent model of galaxy 
evolution. 

In the remainder of the paper we assume that these values 
for (e*,e w ) hold for any galaxy present in the MT. This might 
a poor approximation, but given the persisting ignorance on star 
formation and feedback processes a better performance of more 
complex choices is not guaranteed to produce a more solid result. 
In any case the model success in reproducing key Galactic prop- 
erties provides at least a first and necessary step towards more 
refined treatments. 



3.2 Galactic outflows 

We illustrate the properties of galactic outflows by focusing on 
the case of the progenitors of a 10 13 Mq halo at redshift z = 
(i.e. a small cluster/group). Since the total number of galaxies 
at any intermediate redshift is very large we will deal with mean 
quantities to describe the global properties of the winds. This pre- 
cludes the inspection of one-to-one relationships between galaxy 
and its wind properties, but will make possible to appreciate the 
correlations among several global quantities, as we will see in the 
following. 

Fig.[3]gives an concise overview of the bubbles we have iden- 
tified around galaxies at selected redshifts (z = 10, 6, 3, 1) and 
it can be used to elucidate many physical aspects. Bubbles tend 
to become older towards lower redshifts with a decreasing age 
spread; their ages tend to accumulate close to the Hubble time. 
This behavior reflects the fact that after an early evolutionary 
phase in which bubbles grow by number around relatively low- 
mass objects, the average growth of the halo population mass 
associated with deeper potential wells, prevents the formation of 
new bubbles after z = 3. Hence most of the bubbles seen around 
galaxies at intermediate redshifts (z = 3 — 5), where they are more 
easily observed, contain material expelled by galaxies at a much 
earlier time. At the same time old bubbles merge together to form 
a fewer large ones. The role of merging events is made clear by 
the solid curves in the top panels which show the distribution of 
merging events for halos that have witnessed at least 1 merging 
event. A fraction M/N, indicated by the labels inside the panels 
in Fig. [3] have evolved in isolation until the considered z, i.e. they 
never suffered a merger. Such ratio decreases steadily with time, 
going from 81% at z = 10 to at z = 1. The merging activity 
increases the average ages of bubbles: as discussed above, at high 
z a large number of young, isolated bubbles exist that are later 
turn into aged, large ones. 

The typical bubble size grows from 5 kpc (physical) at 
z = 10, to about 100 kpc at z = 1. Merging is definitely an 
important agent behind this evolution; however, it cannot be the 
only growth factor for bubbles. This can be realized by combining 
the information from the age and size distributions. At the lowest 
redshifts shown in Fig. [3] the bubble age distribution essentially 
collapses onto the cosmic age value; nevertheless, a considerable 
size spread, ranging from 1 to a few hundred kpc, exists. Such 
spread is caused by the corresponding dispersion of galaxy star 
formation rates. Larger galaxies tend to have faster winds, which 
show up as a tail in the PDF of the shock velocity distributions 
(V s > 30 km s _1 ). Hence a minority of the bubbles found at z < 3 
are constituted by old bubbles which are still powered by their 
current galaxy, which is undergoing vigorous star-formation. In- 
stead of creating a new bubble, these massive cluster progenitors 
blow their winds into pre-existing, old bubble structures. 

From the velocity PDF evolution, it is interesting to note 
that at early times bubbles (somewhat contrary to naive expec- 
tations) expand on average at larger velocities and they tend to 
decelerate at lower redshifts. This is an obvious implication of 
the aging trend described above. As larger galaxies dilute their 
energy on the larger volume of pre-existing cavities, they can sus- 
tain the growth of bubbles but cannot re-accelerate them to the 
high velocities at which they were traveling in the past. 

The evolution of the bubble temperature Tj, (Fig. [3j is gov- 
erned by competing effects. Adiabatic expansion and radiative 
cooling (depending on bubble gas temperature and, to a much 
limited extent, gas metallicity) lead to a decrease of Tj,, while 
thcrmalization of the wind energy provide a heat input. Almost 
all the bubbles have temperatures > 10 5 ' 5 K at z = 10, whereas 
by z = 3 this fraction decreases to < 20%. This is mainly due to 
the adiabatic expansion of the bubbles. The expanding outflow is 
powered by the continuous injection of SN energy both in kinetic 
(in practice stored in the shell, where most of the mass is located) 
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Figure 3. Probability Distribution Function (PDF) for different bubble properties (age, bubble radius Rg, shell velocity Vg and internal 
temperature TJ,) at different rcdshifts (left to right) z = 10, 6, 3, 1; error bars show the la dispersion. The solid curves in the top panels 
show the distribution of the number of mcrgings for the halos in the MT. The label M indicates the number of halos that evolved in 
isolation up to that redshift; TV indicates the total number of halos at that z. The dashed lines in the top panels marks the age of the 
universe at that redshift. The dashed lines in the third row panels marks the IGM sound speed at that redshift. 
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Figure 4. Top: Mean kinetic (blue) and thermal (red) energy of 
regions affected by outflows at different rcdshifts. Bottom: Volume 
filling factor of outflows. Shaded areas show la dispersion around 
the mean (solid curve) of fOO merger tree realizations. 



and thermal (stored in the bubble interior) forms. More quanti- 
tatively, from Fig. [4] one can appreciate that the kinetic energy 
is always dominant (by a factor 2-3) with respect to the internal 
one. This kinetic energy will be rapidly converted in disordered 
motions by instabilities, finally resulting in a fully developed tur- 
bulent spectrum. The bubble volume filling factor also increases 
with time, as illustrated in the bottom panel of the same Figure, 
reaching about 10% at z < 2. One has to keep in mind that our 
study is based on the MT of a cluster /small group and therefore 
it might reflect a somewhat biased cosmic regions, although we 
do not expect major differences as we extrapolate these results 
to general IGM. A similar treatment of the out flow evolution has 
been successfully used bv lBertone et ah and ISamui et al.l 

(2008) to in vestigate the r ole of galaxies in enriching the IGM. 
In particular Sa mui et al. use the modified PS formalism 

to follow the evolution of halos and porosity- weighted averages to 
investigate the global influence of winds. For a reasonable choice 
of the model free-parameters, they find that outflows can gen- 
erally escape from the low mass halos (M^ < 10 9 Mq) which 
then dominate the observed enrichment; nevertheless, MW-type 
galaxies ca n create Mpc-s i ze me tal enriched bubbles in their sur- 
roundings. iBertone et al.l l|2005h apply the semi-analytical pre- 
scription for galactic winds to high-resolution N-body simulations 
of field galaxies. Their main conclusion is that galactic outflows 
do not perturb the structure of the Lya forest, galaxies with 



10 9 < M* < 10 10 Mq are the main responsible for IGM chemical 
enrichment at z = 3. 

The bubble volume filling factor is a useful tool when ap- 
plying our models to quantify turbulent effects in the quasar ab- 
sorption line spectra, as we will try to do later on. As a note of 
caution, the filling factor shown in Fig. [4] might be overestimated 
to some degree as we do not attempt to account for possible, 
although likely rare, overlaps between bubbles. 



3.3 Turbulence evolution 

The best way to observationally probe the turbulent content of 
gas affected by outflows is to measure the Doppler parameter, 
b, of absorption lines arising from kinetically perturbed gas. For 
this reason in the following we will use b to quantify turbulence 
evolution. 

In general, one can write b as the quadratic sum of the ther- 
mal and turbulent components, 

& = (<4 + fc?) 1/2 , 



(23) 



vherc 



f 2k B T \ 
\ m a J 



1/2 



is the thermal contribution and 



(24) 



(25) 



is the contribution due to the turbulent pressure defined in Eq.l22l 
Note that the turbulent contribution (contrary to the thermal 
one) does not depend on the mass m a of the atomic mass of 
the absorbing element. This allows to generalize our results to all 
species associated with the cold HI shell gas (e.g. CIV). More- 
over, since the observed Doppler parameter is the sum of the two 
contributions and the thermal broadening is smaller for heavier el- 
ements, the relative contribution of turbulent broadening should 
more important (and hence more easily detectable) for heavier 
elements. 

In Fig. \5\ we show the distribution of Doppler parameters 
in simulated galaxies at different redshifts. The bt distribution 
function at all redshifts shows a similar shape, with a low-energy 
cut-off, followed by a maximum and a steed declining tail towards 
larger values. However, the PDF shift towards lower bt values with 
time: while at z ^ 6, the region affected by the winds can become 
very turbulent, reaching in a few spots bt > 50 km s" 1 , as a result 
of the violently expanding of fresh new bubbles carved by early 
galactic winds, at later times winds blowing in pre-existing bub- 
bles limit become more gentle, producing in decrease of turbulent 
broadening effects. By z « I — 2 the distribution has stabilized on 
median values of the order of 1-2 km s _1 , and even the most tur- 
bulent spots do not exceed 25 km s _1 . The trend above is further 
enhanced by the dissipation of turbulent energy with time, which 
does not allow to store kinetic energy into eddies indefinitely. 

To see this quantitatively, it is useful to estimate the turbu- 
lent dissipation timescale. Under the hypothesis of a fully devel- 
oped, steady-state spectrum the dissipation timescale is well ap- 
proximated by the eddy turn-over timescale, namely t e ~ R s /V s . 
If we further assume that to the early stages of bubble evolution 
the Sedov solution applies, it follows that t e ~ 2.5tb, where i;, is 
the bubble age, in turn shorter than the Hubble time. This im- 
plies that, due to their high velocity and relatively small radii, 
dissipation acts very efficiently at earlier stages, thus decreasing 
the turbulent energy budget of these objects. 

We pause briefly to emphasize an important point. One has 
to keep in mind though (see Fig. [4j that the filling factor of the 
turbulent regions is quite small at high redshift, and therefore it 
must be made clear that the above is not the bt parameter appli- 
cable to the global IGM, but only to tiny and specific (but well 
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Figure 5. Turbulence properties in terms of the Dopplcr parameter bt at four selected redshifts: z = 10,6,3,1 (left to right columns). 
The first row contains the PDF of b t , the others (top to bottom) show the relations between bt and stellar mass, M,, hydrogen column 
density, Njjj, shock radius, R s , respectively. Error bars denote la dispersion in each bin. 
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consequence of the fact that the SFR (which is ultimately pro- 
portional to the energy input rate due to supernovae) is more 
sustained in larger galaxies. However, we see that the onset of a 
turbulent regime is almost unavoidable in regions that are stirred 
by galactic winds. Again, the peak at very low M* is interpreted 
as the imprint of young bubbles in which turbulence has not yet 
suffered from dissipation: we double-checked this hypothesis by 
studying the correlation in different bubble age bins. This anal- 
ysis is reported in Fig. (for a single MT realization) whose 
examination clarifies that the youngest bubbles have the largest 
bt values. Thus, the bursting star formation activity going on in 
small galaxies momentarily reverses the trend of increasing bt as 
a function of stellar mass, clearly indicating that the level of tur- 
bulence depends both on age and mass of the object. The spread 
of the distribution is very limited, suggesting that different mass 
accretion histories do not influence substantially final Doppler 
parameter distribution. 

A relatively easily accessible observational quantity when 
studying the IGM is the neutral hydrogen column density, Nbi> of 
the absorbing shell/system. A reasonable estimate for such quan- 
tity can be obtained from the following formula: 



N,. 



r 

JRs- 



nj](r)dr S3 njjAR s 



Ms 



-AR S 4?r^ij2 
In the latter njj is the mean density in the thin shell, 
the ratio between its mass (M s ) and its volume V s : 



Vs 



4TrRiARs 



(26) 
given by 

(27) 

The neutral hydrogen ionization correction can be obtained 
by solving the appropriate photoionization equ ilibrium bal- 
ance equations given, e.g. by iTheuns et al.l lll998T ) and assum- 
ing the UV backgroun d history of the Universe computed by 
lHaardt fc M adau (1996]). The third row of panels in Fig. |51 shows 
the distribution of bt with N^i- One can note that at redshift 
z ~ 3 the turbulent broadening starts making an impact in HI 
lines with N(HI) ^ 10 14 cm" 2 . 

Finally, the lowermost panel in Fig. [5] shows the cor- 
relation between bt and the distance reached by the out- 
flow. This migh t have important observational implications. 



lAdelberger et al.l d2003h by studying the correlation between 
galaxies and CIV systems with large velocity spreads, showed that 
most of these systems are associated with star-forming galaxies. 
It is conceivable to interpret such result as a direct evidence of 
large-scale outflows around 2 ,$ z Si 3 galaxies. A similar investi- 
gation, as discussed in the next Section, could directly probe the 
level of turbulence in these objects and test our predictions. 



Figure 7. The Doppler parameter bt as a function of the galaxy 
stellar mass for a single realization of the MT, shown at z = 6 
(red circles), z = 3 (green circles), 2 = 1 (blue circles). The filled 
points indicate bubble whose age is less than half of the Hubble 
time. 



defined) portions of it. By plotting the volume-weighted evolu- 
tion of (bt), shown in Fig. [6] we can find support to the previous 
statement: at higher redshift the turbulent pressure contribution 
is less important because of the small size of outflow-affected re- 
gions, whereas in local objects the bubbles grow enough to oc- 
cupy a significant fraction of the volume; however, this occurs 
too late, when dissipation has already started to quench eddies 
on the smallest scales. The trade-off between injection and dissi- 
pation, conspire to produce a maximum in (bt) at z Ri 1. 

Let us turn back to the analysis of the other panels in Fig. [5] 
It is interesting to determine the level of turbulence in bubbles 
produced by galaxies of different stellar masses. In general, bt is 
positively correlated with the stellar mass of a galaxy. This is a 



4 TEST THE MODEL AGAINST DATA. 

A direct comparison of our results with the actual level of tur- 
bulence in the IGM is not possible since this observable has not 
been convincingly derived yet. In the meantime, it is possible to 
perform at least a sanity check for our model as follows. 

In order to study the properties of the IGM, a number of 
high-quality QSO absorption line experiments has been carried 
on during the last decade allowing the build-up of a large sample 
of observed systems for which both TVjji and the Doppler parame- 
ter b have bee n fitted (see an almost complete list of references in 
lMeiksinll2009t) . The observed HI column densities cover the range 
10 12— 22 cm -2 and follow an almost perfect power-law distribu- 
tion dN/dN m oc N~f with = 1.5 - 1.7 . The Doppler parame- 
ters of these systems range over about 10 < b < 100 km s -1 , with 
the vast majority being concentrated between 15 — 60 km/s. The 
observed broadening could be either real, i.e. caused by physical 
mechanisms (e.g. thermal broadening, turbulence and peculiar 
motions), or artificial, due to the blending of close systems. 

Many authors have suggested the use of the lower-cutoff in 
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Figure 8. Left panel: Doppler parameter - HI column density relation observed in the IGM at z = 3 (data points from VLT/UVES 
LP). The red curve represents a fit to the lower cutoff used to derive the thermal contribution. Right: Turbulent Doppler parameter 
distribution at z = 3 as predicted by our model (histogram), also reporting 1-<t dispersion. The red arrows show the observational upper 
limits. 



the Doppler parameter distribution vs. Nm as a reliable proxy 
for the gas temp erature at a g i ven r edshift. Similar approaches 
has been used by ISchave et ail (200Ch to report evolution in the 
inf erred IGM ga s temp erature over the range 2.0 < z < 4.5 and 
by iRicotti et alJ l|2000h to extract the equation of state of the 
IGM from the Doppler parameter distribution of the low-density 
Lya forest. To make progress we follow this procedure. 

In Fig.|8]wc show the absorber properties derived in a sample 
of Lya forests from 22 high-resolution VLT/UVES Large Program 
QSO spectra jD'Odorico et al.ll200sl ; ISaitta et alj|2008h . selected 
at redshift z ~ 3. To avoid the contamination of absorbers affected 
by the UV background flux from the QSO, we exclude from our 
sample absorbers closer than 3000 km/s to the QSO and we also 
do not consider in our analysis Doppler parameters > 100 km s _1 
since they are very likely due to errors in the fitting procedure. 
In the same graph we report our estimated thermal Doppler pa- 
rameter (6th) as a function of iVjji, obtained from a linear fit of 
the minimum b-values, defined by dividing the Nm range in 50 
bins and selecting the smaller b-value in each bin. 

The sanity check then consists in the comparison of our pre- 
dicted mean bt distribution with the analogous average residual 
Doppler parameters from the data after the thermal contribution 
has been subtracted. Given the reasons for contamination listed 
above what we show in Fig. [8] cannot be intended as a direct 
theory-data comparison as the observed broadening still contains 
a contribution from IGM peculiar motions, which are not directly 
related to turbulent motions. Hence, the data points must be con- 
sidered as strict upper-limits which our theoretical Doppler pa- 
rameters should not exceed. Quite encouragingly, this is indeed 
the case: our model is consistent with the experimental upper 
limits in the entire HI column density range. The uncertainties in 
the exact observational determination of bt does not allow a more 
meaningful comparison at this stage. This will become possible in 
the future if Rauch-type experiments using (lensed) quasar pairs 
will be carried out. The increased availability of such systems, en- 
abled by the SDSS, has been already shown bv lKirkman fc Tvtlerl 
(2008), who have discovered 130 close pairs of QSOs in the SDSS 
DR5. 



As an alternative, the presence of the metal lines in the ab- 
sorbers could also be used as a possible way to obtain an indepen- 
dent measurement of the IGM turbulent level. The measurements 
of the Doppler parameters from absorption features correspond- 
ing to two or more elements (typically, CIV and SIV) assumed 
to be co-spatial, allows the separation between thermal and kine- 
matic broadening contributions. As already mentioned in the pre- 
vious Section, this technique could disentangle the broadening 
due to the kinematic from the thermal one which is instead mass 
dependent. To achieve this result high resolution spectroscopy is 
required in order to properly fit the profile of the absorbers along 
the QSO line of sight and to determine the redshift of the lines 
with sufficient precision to make sure that the ions share the same 
absorber environment. We plan to combine these type of obser- 
vations with our model in a forthcoming work. 



5 SUMMARY AND IMPLICATIONS 

Turbulence is an important physical process in essentially all 
astrophysical environments, from planet formation and atmo- 
spheres, to (proto-)stars and galaxies, but particularly on the 
scale relevant to cosmological phenomena it has received a rel- 
atively meager attention. This is partly understood from the dif- 
ficulty to collect reliable data as observations aimed at the study 
of these aspects are very challenging. Even more surprisingly 
though, the well ascertained fact that the IGM is polluted on 
large scales by heavy elements, carried by powerful winds pow- 
ered by galaxies, should raise strong suspicions that turbulence 
must pervade the IGM, albeit at levels yet difficult to understand. 

The present paper aims at making one of the first steps to 
clarify the amount of turbulent energy deposited in the IGM and 
then dissipated through a cascade to smaller scales by galactic 
winds and to characterize its properties. To this aim we have de- 
veloped a simple but physically accurate model based on standard 
ACDM hierarchical structure formation model to identify halos 
hosting galaxies progenitor of a small cluster /group of total mass 
M = 10 13 Mq at z = 0. In addition, the required star formation 
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history in all their progenitors via a merger-tree technique (care- 
fully calibrated to reproduce the Milky Way) is also followed. As 
these galaxies progress along their history, their collective super- 
nova explosion vent gas, heavy elements and, most importantly 
for the present study, energy into the surrounding IGM. The evo- 
lution of these hot intergalactic bubbles is followed by including 
all the relevant physics, enabling us to determine their thermo- 
dynamic properties at best. We couple then this general galaxy- 
IGM interplay scenario (which can be used and adapted to a 
number of ancillary applications as e.g. metal enrichment and SZ 
effect) to a treatment of turbulence evolution using a powerful 
approach based on the s pectral transfer function developed by 
iNorman &; Ferraral l|l99(Jl . In brief, the major advantage of this 
method is a relatively straightforward derivation of the turbulent 
IGM turbulent pressure evolution (see Eg. \22l . 

The main results can be summarized as follows. At z 3 the 
majority of the bubbles around galaxies are old (ages > 1 Gyr), 
i.e. they contain metals expelled by their progenitors at earlier 
times, and have a size distribution at that redshift in the range 
10-100 (physical) kpc. The velocities reach up to 100 km s _1 for 
larger galaxies but a considerable fraction of the bubbles have 
already become subsonic and the shock decayed into a sound 
wave; temperatures in the rarefied bubble cavities are in the range 
5.4 < log(T/K) < 6.5. The bubble volume filling factor increases 
with time reaching about 10% at z < 2. The energy deposited 
by these expanding shocks in the IGM is predominantly in ki- 
netic form (mean energy density of 1 /xeV cm -3 , about 2-3 X the 
thermal one), which is rapidly converted in disordered motions 
by instabilities, finally resulting in a fully developed turbulent 
spectrum. The corresponding mean turbulent Doppler parame- 
ter, bt, peaks at z 1 at about 1.5 km s — 1 with maximum 
bt = 25 km s — 1 . More informative though is the bt distribu- 
tion associated with individual bubbles and the relations with 
the galaxy stellar mass and bubble shell column density and ra- 
dius (shown in Fig. 5). The shape of the bt distribution does not 
significantly evolve with redshift but undergoes a continuous shift 
towards lower bt values with time, as a result of bubble aging. We 
find also a clear trend of decreasing bt with JVm and a more com- 
plex dependence on R s resulting from the age spread of the bub- 
bles. We have finally tried to compare our results with available 
data on the Doppler parameter deduced from QSO absorption 
line studies of the hya forest, but we are facing the difficulty of 
accurately removing the contamination due to the thermal com- 
ponent and the more troublesome peculiar contribution. Alterna- 
tively one can extend our study to make predictions for quasar 
pairs systems and/or using observations involving different metal 
species. We are already working along these directions, which in 
general require additional information of the position of galaxies 
and therefore the implementation of our turbulent model into a 
full-scale hydrodynamic simulations. 

If turbulent is indeed present in the intergalactic medium, as 
our study seems to suggest, it might have relevant implications for 
a number of IGM studies. First of all the long-standing problem of 
a discrepancy between the observed Doppler parameter observed 
and those deduced from numerical simulations of the hya forest 
mig ht be alleviated by a proper inclusion of turbulence broaden- 
ing. lOppcnhcim er fe Pavel Il2009h convincingly demonstrated that 
only a model where explicitly add sub-resolution turbulence, and 
the required turbulence increases with O VI absorber strength, is 
able to match all the observations, giving an hint of this correla- 
tion between turbulence and enrichment that we are planning to 
explore in a forthcoming paper. Secondly, the turbulent energy 
injection might affect the matter power-spectrum determination 
from experiments us ing the Lya forest as a probe, particularly on 
the smallest scales ((McDonald ct al. 2006). Third, if turbulence 
is largely present in the IGM it might result in a detectable sig- 
nature in the scintillat ion of distant quasars ( Fer rara fc Perrial 
l200ll ; lLazio et aHl2004T) . The IGM hosts the necessary condi- 



tions for scattering, is largely ionized and permeated by shocks 
from large scale formation or galactic feedback. Moreover, the 
long path lengths through the IG M can compensate for the lower 
density with respect to the ISM. lLazio et al.l j2008f ) have consid- 
ered intergalactic scattering , expanding on previo us treatments 
from intra-cluster media of lHall fe Sciamal l ll979h . to scintillat- 
ing AGNs, their findings are broadly consistent with the scenario 
of a highly turbulent intergalactic medium. This possibility was 
explored so far only as an Ansatz, to which our study now pro- 
vides a more solid support. Fourth, turbulence might be carefully 
modeled in order to interpret the results for experiments aim- 
ing at directly measuring the change of the expansion rate of 
the Un iverse with time and the variability of f undamental con- 
stants llCorasaniti et all |2007| ; iLiske et al.ll2008h , as for example 
the p roposed CODEX-ESPRESSO experiment llCristiani et al] 
120071) . Last but not least, the relatively mild level of turbulence we 
find at z fs 3 is consistent so far with the upper limits coming from 
observations. Moreover, it indicates that the IGM at those epochs 
had already dissipated a large fraction of the kinetic energy de- 
posited at early times (2 > 5) by the galaxies which initiated the 
metal enrichment and, possibly, the reionization process. All this 
well agrees with the idea that the metals we detect at moderate 
redshifts were mostly produced during a " pre-enrichment" phase 
of the IGM, by a population of small galaxies ancestors of the 
ones we study routinely in detail, the Lyman Break Galaxies. 
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